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CONVECTIVE HEAT TRANSFER FOR SHIP PROPULSION 
M. A. Habib* and D. M. McEligot 


Aerospace and Mechanical Engineering Department 
University of Arizona, Tucson, Ariz. 85721 USA 


Abstract 


Calculations of the flow and heat transfer of a swirling turbulent 
flow after a sudden pipe expansion are reported. The calculations were 
obtained by the numerical solution of the time-averaged forms of the 
continuity, momentum and thermal energy equations together with equations 
for the turbulence kinetic energy and its rate of dissipation. For a 
sudden expansion without swirl, the predictions produced satisfactory 
agreement with available data for the Nusselt number. In the swirling 
case, there were no heat transfer measurements available for comparison. 
However, existing measurements of a swirl flow at an expansion without 
heat transfer were used to test the validity of the flow field calculations. 

The predicted effects of varying the swirl number from zero to 1.0 
on the heat transfer behavior are presented for a range of Prandtl numbers 
from 0.7 to 10. The expected effects of the swirl on the velocity and 
temperature fields are also reported. The results predict, for example, 
that the maximum Nusselt number is increased as the swirl number increases 
and its position moves towards the inlet section. They also show an in- 


crease in the Nusselt number as the Prandtl] number increases. 


*Presently at the Mechanical Engineering Department, Cairo University, 
Cairo, Egypt. - 
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NOMENCLATURE 
Symbol Definition 
Cc. Constant in turbulence model, for correction for stream- 


line curvature, equation (11) 


D Diameter of pipe downstream of the expansion = 2k, 
d Diameter of the pipe upstream of the expansion = 2r 
G Generation of turbulence kinetic energy K 
H Step height = (D = d)/2 
H Enthalpy 
K Turbulence kinetic energy = 5 or 
Nu Nusselt number 
Nue Nusselt number of fully developed pipe flow 
Pry) Molecular Prandtl number 
Pry. Turbulent Prandtl number 
q, Wall heat flux 
R Radius of downstream pipe 
r Radial coordinate 
Re Reynolds number; Regs based on diameter upstream; Rep, 
based on diameter downstream 
Re Flux Richardson number 
$ Swirl number, sf ir ar/r_f*o Trdr 
Ce) @) 
Ty Bulk temperature, {UTrdr 
sUrdr 
T Inlet temperature | 
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| W Tangential mean velocity 
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i € Rate of dissipation of turbulence kinetic energy 
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CHAPTER I 


INTRODUCTION 


1.1 Background 


Current Naval propulsion plants are powered by variations of the 
Rankine cycle (steam) or the open gas turbine cycle (air and combustion 
products) plus some diesel engines in small ships. Alternative power 
systems suggested include the closed gas turbine cycle and cycles involving 
dissociation of the working fluid in either a Rankine or a gas cycle. 

These latter two are believed to offer the potential of substantial im- 
provement in the power-to-weight ratio of the propulsion plant. The present 
and proposed studies consider basic problems in convective heat transfer 

and flow friction that are important in all of the above. 

Convective heat transfer provides the dominant thermal resistance in 
several components of conventional steam power plants as well as in all 
heat transfer components in gaseous cycles. For example, heat transfer 
from the condensing steam to the cooling water in the main condenser of 
a Naval ship is via the condensation heat transfer coefficient which is 
high, then conduction conductance through thin tubing of high conductivity, 
and finally via the moderate heat transfer coefficient of the cooling 
water flowing through smooth circular tubes. Thus, the overall thermal 
resistance is dominated by the convective thermal resistance of the cool- 


ing water inside the tubes. One can expect significant reductions in tube 


length and, therefore, size and weight of the condenser and overall plant 
if the convective heat transfer coefficient of the cooling water is in- 
creased appreciably. 

A common geometry recurring in compact Naval propulsion plants is 
a change in duct size in the primary fluid loop. As a consequence of the 
upstream plumbing, the fluid is often swirling about the axis in the 
piping. Heat losses from the primary fluid and thermal stresses in the 
component depend on the convective heat transfer from the fluid to the 
component as it undergoes this geometrical transition. The idealized 
problem is a study of heat transfer in a sudden expansion with swirl flow, 
the subject of the present report. The same situation also occurs at the 
entrance to some heat exchanger tubes. 
The flow downstream of a sudden expansion consists primarily of a 


central jet, with or without rotation due to swirl, surrounded by a re- 


circulating donut-shaped vortex induced by the separation as the jet 
leaves the central tube. This jet gradually spreads as the flow proceeds 
downstream until it reattaches to the downstream wall and fills the larger 
tube, Eventually a fully developed velocity profile can be expected to be 
established if the tube length is sufficiently long. The centrifugal 
body forces induced by swirl modify the recirculating vortex, the develop- 
ment of the jet and the reattachment location. 

Significant enhancement of heat transfer rates generally occurs 
when fluids are transported through sudden enlargements in pipes (Fig. 1) 


or duct cross-sections or when they are subjected to angular tangential 


r 
Adiabatic wall ~ 


Figure 1. Idealized qeometry. 


momentum due to twisted tapes or the application of vane-generated 
swirl. In both cases, the levels of turbulence kinetic energy are 
increased and the heat transfer rate generally increases as well. 

High shearing rates created by both separation and swirling effects 

are associated with high generation rates for the kinetic energy 

of turbulence; in addition the rate of dissipation of turbulence kinetic 
energy decreases due to an increase in the length scales. The two 
effects lead to elevated magnitudes of the kinetic energy of turbulence. 
All these effects are expected to reduce the thickness of the viscous 
sublayer through which heat must pass largely by molecular diffusion 


and so will increase the heat transfer rate. 


1.2 Previous work 


The problem of heat transfer to separated and reattached sub- 
sonic turbulent flows downstream of abrupt pipe expansions has veen 
studied experimentally by a number of investigators such as Ede, Hislop 
and Morris [1956], Krall and Sparrow [1966] and Zemanick and Dougall 
[1970]. Reviews of the heat transfer literature for separated and 
reattached flows are reported by Eckert et al. [1971] and by Fletcher, 
Briggs and Page [1974]. The experimental data of Krall and Sparrow 
(1965] for water, and of Zemanick and Dougal [1970] for air are for 
typical geometries representing axisymmetric abrupt enlargement of 
tubes. Both teams found a maximum in the local heat transfer coeffici- 
ent at the assumed point of flow reattachment. Adequate descriptions 


of wall heat transfer at solid boundaries in regions of separated and 


BS ec ce a ss ae ee 


reattached flows have been limited by the lack of a complete solution 
for the fluid flow field. Spalding [1967] presented a power law re- 
lation between the Stanton and the Reynolds numbers for heat transfer 
from a wall adjacent to a region of separated flow; he employed a one-~ 
dimensional model of the flow near a wall. Recently, Chieng and Launder 
{1980}, in calculations of the turbulent heat transfer downstream of 

an abrupt pipe expansion, found that the Nusselt number is overpredicted 
by 20% in flows with high Reynolds numbers and that its maximum value 
occurs about one step height upstream of the calculated reattachment 
point, relative to the data of Zemanick and Dougall. 

For swirling flow, a number of studies have been done to investi- 
gate the effect of swirl on heat transfer parameters. Experiments 
performed by Thomas [Eckert et al., 1971], with swirling air flow 
in a tube indicated that, for turbulent flow and Reynolds numbers less 
than 10°, the swirl flow transfers heat more efficiently than a non- 
swirling flow. Zaherdadeh and Jagadish [1975] concluded that an aug- 


mentation of up to 80% was obtained, with a constant wall temperature 


boundary condition, with the use of tangential vane type swirl generators. 


1.3 Present work 


Calculations of the effects of separation and swirl on heat transfer 


parameters have rarely been reported and the combined effects of both 
apparently have not been studied. Thus, we attempt the prediction of 


such cases in contrast to the "postdiction" usually employed in numerical 


OE 


studies, In the present paper, a calculation method is provided and 


. . ae — — , = EA — ss eS esti ant 


is applied for the prediction of flow and temperature fields in the 
swirling, separated and reattached flow in the vicinity of a sudden 
expansion of a pipe. 
An idealized version of the geometry is shown as Figure 1. The 
expansion ratio is 2 and the Reynolds number upstream of the expansion 
is taken as 5 x 10° for the calculations. The effects of swirl on the 
heat transfer parameters are studied in the range of swirl number 
from 0.0 to 1.0 and of Prandtl number from 0.7 to 10. For confidence 
testing, the calculations are compared with the measurements of Zemanick 
and Dougall [1970} for zero swirl number and those with swirl, but without 
heat transfer, are compared to the data of Beltagui and MacCallum [1976]. 
The governing equations and boundary conditions are described 
in the following section which is ended by a description of the solution 


procedure. The results are then presented and discussed. The paper 


ends with a summary of conclusions. 
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CHAPTER 2 | 


ANALYSIS ‘ 


Equations representing conservation of mass, energy, turbulence 
kinetic energy and its dissipation rate plus the momentum equations 
were solved in axisymmetric coordinates with an extension of the 
TEACH computer program developed originally by Gosman and Pun [1974]. 
In the present study, the version of Habib and Whitelaw [1980], which 
treats double coaxial jet flows, was modified to give predictions of 
swirl flow downstream of an abrupt pipe expansion and to yield the 
temperature field and heat transfer parameters by adding the solution 
of the thermal energy equation. 


The equations modeled and the turbulence model are as follows: 


continuity, ~_ 
U: 
2 J = Oo (1) 
dX; 
momentum, 
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thermal energy, 


(3) 


\~ 
ae 
“od 

C 

Lr 
Na 

1] 

ha 
pees 

E 

a 

hes 

i 
YS 


turbulence kinetic energy, 


ao ale eee 


and dissipation of turbulence kinetic energy 


2 (plel= 3 (SS) +ce eS (5) 


ax.\ Se ay; Mk 2k 
cupk~ 
with M ece = a + - nf (6) 
G=-e GG, a , k= @u,/2 
j 


pu; = ere a0: 4 4) - seks, 


where 555 =0 if i#j and Fe = 1 if i= j. (7) 
Fluid properties were idealized as constant except in some cases with 
air where the density was allowed to vary due to pressure variation. 

In addition to the set of equations listed above, the following 
inlet and boundary conditions are considered. In the case of non- 
swirling flows, data for fully developed pipe flow were taken as initial 
values for mean velocity and turbulence kinetic energy. The inlet 


dissipation rate was then determined from its model, 


3/2 
e=c Sy (8) 
u 


where 2 is the length scale that characterizes the energy containing 
eddies, taken as 2% = 0.03 T° 
For swirl, the initial axial and tangential mean velocities were 


interpolated from the data of Habib [1980]. The turbulence kinetic 


Mme sewers. 


hn 


energy was taken as a function of the total mean kinetic energy and the 
dissipation rate was determined via equation (8) with 2% assumed to be 
given as ro (0.03 + 0.078). 
All axial gradients, 3/9x, were presumed zero at the exit plane 
of the tube well downstream. Axial symmetry was specified so V = 0 and the 
radial gradients of all other quantities were zero at the centerline. No 
slip and impermeable wall conditions were specified on all solid surfaces. 
Wall functions [Launder and Spalding, 1972] corresponding to the 


dependent variables were taken as 


un. y* 
epee a) 


+ 1 + + 1/4 ,,1/2 
where U ==—2&nE and =C K 
Re np ee YO KS ¥p/Ms 


+ 
u,(H_ - H_) y 
and oe e (10) 
q, + 
y U_ + Pe 
P P 


(Fe \ Pry, 1/4 "1 
where P. = C. {—" - 1} 
f f \Pr,, ] Pry, 


The value of the kinetic energy of turbulence near the wall, Ky is 
calculated from the transport equation for K with the flux of energy 


to the solid wall set to zero. The corresponding value of « was cal- 


culated from equation (8) with 2 = ce ae 


*) 
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The constants used in the above equations are in accord with those 
of Khalil and Whitelaw [1974] and Pope and Whitelaw [1976]. They are 
summarized as follows: 


Ge = 0.09, o, = 1.0, o. = 1.22, 0. 


K 
c, = 1.9, C, = 9.24, « = 0.42, E = 9.8 


Modification of Cc, to account for streamline curvature is described in 
Section 3.2 later. 

The solution procedure described by Patankar and Spalding [1972] 
was used to solve the above equations. From trial values of the pressure 
distribution, the momentum equations can be solved to obtain an esti- 
mate of the velocity field at the first iteration. However, there is 
no guarantee that the resultant velocity field will satisfy the continu~ 
ity equation, therefore, after each calculation over the solution 
domain, adjustments are made to the pressure and velocity fields to 
satisfy continuity. The procedure is repeated until convergence is 
obtained. The energy equation is then solved and the temperature field 
is presented. 

Most of the computations were done with a grid of 25 x 18 nodes 
with more points concentrated near the walls and regions of separation 
and high velocity gradients. The computer time (per iteration per 
grid node) using the CDC CYBER 175 was 9.5 x o-9 Sec to solve the velocity 
field and 1.2 x 1074 Sec to obtain the temperature field. (These time 
estimates typically corresponded to $24/run and $6/run, respectively, 
at University rates). The storage requirements of the program were 


37,000 octal words. 


ll 


A grid independence test was conducted during the comparisons to 


the data of Zemanick and Dougall [1970]. It showed that an increase 


9 or oat eta nasite: nieteamealenilen aetna theta debaeriatiilallay 


in the number of nodes from 20 x 17 to 27 x 21 affects the predicted 
Nusselt number only in the region of its maximum (see Figure 2) and 

then less than three percent. The effect was negligible in the down- 
stream region. For most of the predictions presented in the next section 


a grid of 25 x 18 nodes was employed so they are expected to have better 


numerical accuracy than the former case. 
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© Zemanick and Dougall [1970 | 


Chieng and Launder [1980] 


Re= 47,640 
d/D= 0.43 


Figure 2. Comparison for heat transfer at a sudden expansion without swirl. 
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CHAPTER 3 


PREDICTIONS 


In this section results from the computer program are first compared 
to the available data to test its validity. Then, the predictions of 
the separated flow with both swirl and heat transfer are presented. 
3.1 Zemanick and Dougall - heat transfer without swirl 

The results are first compared with the experiment of Zemanick 
and Dougall [1970] for zero swirl. Since inlet velocity profiles were 
not presented, the initial conditions were taken as those of fully 
developed pipe flow in the smaller tube. In the run chosen for compari- 
son, the expansion ratio was 0.43 and the Reynolds number calculated 
downstream from the expansion was 47,800. In the experiment resistive 
heating was employed so the appropriate thermal boundary condition would 
be a specified wall heat flux. To ‘eseveese with the procedure of the 
present program, the wall temperature distribution which was measured 
was entered as the thermal boundary condition and the wall heat flux 
was calculated via equation (10). 

As noted by Zemanick and Dougall [1970] there were some difficulties 
with their experiment so it is not ideal for a close comparison, but 
it is one of the few with data readily available. Axial conduction 
in the tube was significant and the thickness of the copper electrode 
at the expansion extended about a step height in the downstream direc- 


tion. In addition it may be noticed that their correlation of local 
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Nusselt numbers in the fully developed region downstream is about 15% 
higher than commonly accepted values [McEligot, Magee and Leppert, 
1965]. Unfortunately, quantitative estimates of the experimental 
uncertainties [Kline and McClintock, 1953] were not presented. If 
one assumes an uncertainty of 1°C in measuring the wall temperature, 
typical of premium grade thermocouple wire, it can be shown that this 
effect alone would lead to an uncertainty of about twenty percent in 
the Nusselt number near the entry and approximately five percent down- 
stream; these estimates are shown by brackets on Figure 2. Thus, it is 
best to restrict the comparison to trends rather than absolute magnitudes. 
Figure 2 shows a comparison between the calculated and measured 
distributions of Nusselt numbers; the curves represent the numerical 
predictions. It is shown that the agreement is generally good in the 
upstream and downstream locations with the present calculations under- 
predicting the maximum Nusselt number by about 20%. The location of 
. the maximum Nusselt number is found to be one step height upstream 
from the corresponding peak in the experiments. By examination of 
the predicted velocity distributions it may also be found that the 
calculated maximum Nusselt number occurs two step heights upstream 
from the calculated reattachment point. This direction was also indi-~ 
eated by Chieng and Launder [1980]; however, they indicate the 
difference in location to be one step height only. 
The discrepancy between the calculations and the experiment may 


| be attributed partly to the numerical accuracy in the calculation 
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procedure. Increasing the number of grid points to reduce effects of 
the numerical diffusion leads to a slight increase in Nusselt number 
in the entry region as shown on Figure 2. However, as noted earlier, 
the increase in the number of nodes from 20 x 17 to 27 x 21 gives a 
change of less than 3% which is essentially negligible for the present 


purposes, 


3.2, Beltagui and MacCallum - adiabatic flow with swirl 


For swirl, data with heat transfer are not available to compare 
to predictions; however, the adiabatic measurements of Beltagui and 
MacCallum [1976] were considered to test the flow field calculations 
and to provide guidance. It is expected that when the model provides 
predictions which are satisfactory for the flow field calculations, then 
the heat transfer calculations which depend on them will consequently 
be reasonable too. 
Predictions were calculated for the conditions of the experiment 
by Beltagui and MacCallum [1976] with Re = 2 x 10°, S = 0.67 and 
d/D = 0.413 as expansion ratio. Their measurements do not include 
inlet velocity profiles; therefore initial values for the present cal- 
culations were interpolated to S = 0.67 from the data of Habib [1980]. 
Results based on the standard K-e model of Section 2 are shown 
as curves with centerline markings (dash-dot) in Figure 3(a) for axial 
mean velocity profiles and Figure 3(b) for circumferential mean velocity 


distributions. Dashed lines represent the experimental observations, 


It should be noticed that the ordinate scale changes as effects decay 


downstream. 
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One effect of this relatively high swirl is that the central jet 
rapidly takes the form of a circular sheet which moves radially outward 
and attaches to the wall of the outer tube in a short distance. Conse- 
quently, the surrounding vortex becomes significantly smaller in size 
than without swirl and is seen as a small reversed flow region at only 
the first axial location. Along the centerline of the larger tube 
a reversed flow also appears showing the existance of another recir- 
culating vortex of size comparable to the diameter of the smaller tube. 
The qualitative trends are in agreement with the data. The axial vel- 
ocity profiles agree well except near the centerline where the vortex 
is predicted to be shorter than measured (e.g., compare centerline 
velocities at X/D = 2.63). Also the tangential velocity in the sential 
region is predicted to decay more quickly than observed. 

The discrepancies are partially explained by Khalil and Whitelaw 
(1974] and Habib and Whitelaw [1980], who showed that the standard 
K-e model often yields underpredicted lengths for the central recircu- 
lation zones of swirling flows. The work of Pope and Whitelaw [1976] 
for flows behind a baffle and of Habib and Whitelaw [1981] for confined 
flows with a sudden expansion indicates that the dissipation equation 
is the weakest of the governing equations since its modelling is only 
based on assumptions. As shown by Bradshaw [1973], body forces due to 
streamline curvature can have a strong effect on turbulence and the 


equations of turbulence models should account for these effects. Re- 


cently, Rodi [1979] and Morse [1980] found that it is essential to 


| 
| 
i 
| 
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correct one of the empirical constants in the dissipation equation as a 
function of Richardson numbers to achieve agreement with experiments. 


To improve the present calculations, we replaced C the empirical 


nes 
constant in the dissipation equation, by c,(1 + CR Ry) where Re is the 


flux Richardson number [Bradshaw, 1973] 


(11) 


The use of Cc. = 0.2 was found to increase the length of the recircu- 
lation zone from about 1.7D to 3.6D. This value was chosen to pro- 
vide the best agreement with the measured size of the central recircu- 
lation zone. 

Predictions with CR = 0.2 are presen..d as solid curves in Figure 3. 
The axial velocity profiles show considerably better agreement with the 
data in the central recirculating region, as would be expected from the 
manner of adjustment, and slightly better agreement near the wall. 
The predicted circumferential velocity distribution is also improved 
near the wall, the region with the greatest effect on the heat transfer, 
but as the centerline is approached the values appear to be under- 
predicted. Comparison of the predicted velocities with Ce = Q and 


= 0.2 shows that modifying C, with this functional dependence on 


CR 1 


Re tends to decrease the rate of diffusion in the downstream sections 


and to suppress the development of the jet. 
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3.3 Heat transfer with swirl at a sudden expansion 


Calculations of heat transfer to a swirling flow were made after 
an abrupt enlargement of D/d = 2 for an inlet Reynolds number (Re) 
of 50,000 and with constant wall temperature imposed after the expansion. 
It was idealized that there was no heat transfer through the annular step 
between the inner and outer tubes. For flows at the different swirl 
numbers studied the inlet distributions of velocities and kinetic energy 
of turbulence were interpolated from the data of Habib [1980]. The 
effects of varying swirl number and Prandtl number were examined sep- 
arately. 

Effects of swirl were studied with the Prandtl number taken as 
0.7 as for common gases. Predictions are shown in Figures 4 through 7. 
Figure 4, presenting the development of velocity and temperature dis- 
tributions, provides information which helps understanding the later 
figures. On these figures the solid curves represent the reference 
case; zero swirl. 

In Figure 4 the velocity profiles show that at S = 0.1 the flow 
is only slightly modified, at S = 0.4 there are significant differences 
and for S = 1 the general pattern of the flow changes. These runs 
might be considered as showing low, medium and high swirl, respectively. 
As the swirl (or angular momentum) increases, the reattachment point 
of the separated region in the corner moves successively closer to 


the entrance and the central jet spreads more rapidly towards the wall. 


After a short distance the velocity profile for S = 0.4 resembles a 
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Figure 4, Develonoment of mean velocity and temnerature profiles for 
different swirl numbers. 
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plug flow. The velocity near the heated wall increases with S in this 
entry region. 

At S = 1 the flow pattern resembles that of Beltagui and MacCallum 
{1976] shown earlier. For this high swirl number the flow exhibits 
a long central recirculation zone while the outer recirculation zone is de- 
creased to 0.13H. Flow development is delayed as the central jet be- 
comes a wall jet with the reversed flow at the centerline persisting to 
beyond X/H = 17. These features are also demonstrated by Habib [1980]. 

The corresponding variation of the non-dimensional wall shear 
stress or friction factor is shown on Figure 5 (note the change of 
scale for S = 1). Negative values occur due to the reversed flow near 
the wall. The value of zero indicates the location of the reattachment 
point or position where the dividing mean streamline meets the wall. 
As S increases this point progresses upstream as noted from examination 
of the velocity distribution; the points of maximum shear stress like- 
wise move upstream. The magnitudes of the extremes also increase with 
swirl; an effect of the centrifugal forces apparently is to drive the 
higher velocity central flow towards the wall, thus steepening the 
velocity gradient. For fully developed, turbulent flow in a tube at 
Re) = 25,000 the friction factor is expected to be about 0.006 [Kays, 
1966] so all the flows must go considerably further downstream before 
becoming fully developed. 


The temperature profiles in Figure 4 show approximately isothermal 


zones across the recirculating regions, thereby demonstrating high transnort 
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Figure 5. Development of wall shear stress for swirl flow. 
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rates due to the turbulent mixing in these regions. For the range 


S < 0.4 there is only a slight change in temperature distributions as 
the swirl is increased. 

The interesting case is for high swirl flow. Within 5 step heights 
(or 1.25 D) the central recirculating flow is approximately isothermal 
By X/H = 1 this reversed flow has carried thermal energy forward from further 
downstream and, thereby, has raised the centerline temperature sub- 
stantially. The small recirculating region near the wall has done the 
same for that region, giving the appearance of a cooled wake at an 
intermediate radius. However, before X/H = 5 the central region has 
expanded and has eliminated the wake appearance in the temperature 
distribution. 


Predicted local Nusselt numbers are shown in Figure 6 normalized by 


8 P 0.4 


Nu = 0.020 Re?’ r 


fd 


a correlation for fully developed turbulent flow with a constant wall 
temperature [Kays, 1966]. As for the wall shear stress, increasing 
swirl increases the maximum value and moves it upstream towards the 
entry. The higher velocity gradient corresponds to flow at a higher 
Reynolds number so the viscous layer, which dominates the thermal 
resistance, is thinner and the Nusselt number is increased. Comparison 
to Figure 5 shows that for each swirl number the maximum Nusselt number 


occurs near the point of reattachment, which was shown earlier to move 


upstream as swirl increases, 


Re4= 50,000 
D/d= 2, Pr=O7 


Figure 6, Effect of swirl on development of Jocal ‘usselt numbers. 
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For the high swirl case, the increase in Nae ae due to the combined 
effects of separation at the annular step and swirl is four times while 
for separation alone the improvement is 2.5, relative to fully developed 
flow. Thus, the major increase is due to the separation. Further 
downstream, the low and moderate swirl cases show a slight decrease 
in Nusselt number as swirl number increases. This result may be the 
consequence of swirl reducing the length of the outer separated zone 
so that the development of the normal tube flow has progressed further 
at a given location, i.e., for given X, X - A og sbeacnnene increases 
with S. For S = 1 the strong wall jet phenomenon continues to dominate 
the velocity profile and, apparently, maintains a slightly higher 
Nusselt number than for zero swirl. 

For applications such as heat exchangers, the designer is interested 
in the resulting bulk temperature variation of the fluid as it is heated 
(or cooled). Figure 7 presents this information and it is seen that 
an effect of swirl is to increase the bulk temperature and its rate 
of increase slightly. After viewing the temperature profiles, it may 
seem surprising that the increase is not greater for high swirl. The 
difficulty is that, though the temperature is approximately uniform, the 
velocity changes sign in the large reversed flow in the central region. 
| In calculating the bulk temperature via its definition (/UTrdr//Urdr) 
the enthalpy flow rate there is subtracted from the enthalpy flow rate 
in the positive direction giving a net value which is smaller than 
| expected. Therefore, while local heat transfer rates are increased 

substantially by swirl the overall improvement is not proportionally 


[ as large. 


Re, = 50, 000 
D/d=2, Pr=0.7 


Figure 7. Effect of swirl on overall heating. 
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| Predictions of the effects of Prandtl number on heat transfer were 
calculated for S = 0.7, again at Re, = 50,000 and D/d = 2, Examination 
of the expected velocity profiles in Figure 8 shows this situation to 
be a high swirl case with a substantial recirculating vortex in the 
central region. In comparison to the profile for S = 1 at X/H = 1 
in Figure 4, one finds the same shape except the maximum does not fall 
as close to the wall in this case with lower swirl. Lower circumfer- 
ential momentum at the entry causes the main streamlines of the annular 
jet to have less rapid divergence from the axis. 

Temperature profiles were calculated for Pr = 0.7, 3 and 10, corres- 
ponding to common gases and liquid water at about 60°C and 10°C, 
respectively. Results are also shown on Figure 8. The general behavior 
is in agreement with the earlier predictions for S = l. 

At each axial location the temperature profiles have the same 
general shape for all Prandtl numbers. Only the magnitudes change 
significantly. Thus, the total thermal energy transferred to the flow 


is greater as the Prandtl number decreases. The similarity of the 


shapes may be taken as an indication that convection dominates the 
thermal energy transfer processes in these regions despite the in- 
creased thermal conductivity implied by lower Prandtl numbers. The | 
levels of the curves then depend on the heat fluxes at the wall where 
the thermal conductivity does dominate the thermal resistance. Thus, 
accurate prediction of heat transfer parameters requires accurate 
knowledge of the behavior of the viscous layer in this flow as well as 


in simpler ones. 


T-T, 


in 
TwoTin 


Figure 8. 


Effect of Prandt] number on temperature distributions. 
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The axial distribution of local Nusselt numbers is presented as 
Figure 9; the curve for Pr = 0.7 is consistent with the trends of 
Figure 6 earlier. Again, as with the temperature profiles, the shape 
of the Nusselt number distribution is primarily dependent on the flow 
parameters-~Re, S and D/d--and is relatively independent of the Prandti 
number. However, the level or magnitude of the results increases as 
the Prandtl number increases so, as a general conclusion, one may 


note that the local Nusselt numbers increase with the Prandtl number 


as in most flows. 
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D/d=e, S$=0.7 
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Figure 9. Effect of Prandtl] number on local ‘Nusselt number, 


CHAPTER 4 


CONCLUSIONS 


The following important conclusions can be extracted from the 

preceding results: 

1) For heat transfer without swirl, generally good agreement 
is obtained between measured and calculated Nusselt numbers. 

2) For an adiabatic flow with swirl, predictions could be ob- 
tained in agreement with available data by modifying the 
turbulence model to account for effects of streamline cur- 
vature. 

3) For swirl flow with heat transfer, (a) the maximum local 
Nusselt number is predicted to increase and its location 
moves towards the entry as the swirl number increases, and 
(b) the local Nusselt number increases as the Prandtl number 


increases, 


Comments concerning parameters used in Block data of Teach T.CDC 
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APPENDIX. COMPUTER PROGRAM 


This parameter indicates continuation of calculations 


from previous output saved on magnetic tape 
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Does not read or write any output 
Writes data on tape 8 
Reads from Tape 7 and writes on Tape 8 
Parameter controls heat transfer calculations 
There is a heat transfer calculation (the thermal 
energy equation is solved) 


No heat transfer calculations 


Parameter used to restrict the program to calculations for a 
constant property fluid 
Flow field calculations only 


Energy equation is solved alone provided that the flow 


is known from a previous set of iterations. In this 
case (i.e., ICP = 1) IHT should be set to 1 also. 


The program may consider changes in density when the 
fluid is air. For this case set [HT = 1 and ICP = 0. 


Defines North wall (i.e., larger tube) as adiabatic 
There is heat transfer through the North wall 
There is no heat transfer through the North wall 


Sidewall (end wall) boundary situation 


There is heat transfer through the side wall. 
(This option was not tested.) 


There is no heat transfer through the side wall 


*Questions concerning the details of the program should be directed to 
Dr. M. A. Habib, Power Section, Mechanical Engineering Department, Cairo 


University, Cairo, Egypt. 
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0.2 


WW 


MAXIT 


INDPRI 


PO 


CP1 


GAS CON 


WFLUID 


COND 


NOTAIR 


CPFL 


VISFL 


CONDFL 


ICHF 


A-2 
Constant C, used with equation (11) for streamline curva- 
ture correction. (It is only useful in the swirl cases.) 
There is no correction 
There is a correction (with CR = 0.2) 


Maximum iterations allowed. (It will write the results 
if IGET>2) 


Defines the number of iterations after which the program 
shows a sample of the calculations 


Atmospheric pressure, N/m? 

Specific heat, J/Kg C° 

Universal gas constant = 8314 J/Kgmole C° 
Molecular weight, Kg/Kgmole 


Thermal conductivity, Watt/m C° 


The fluid is air or a comparable gas 


The fluid is not air and the properties for the fluid 
must be inserted. These are identified as CPFL, 
CONDFL, VISFL, DENFL) 


Laminar and turbulent Prandtl numbers 


The properties of the fluid to be used if it is not con- 
sidered to be air. (NOTAIR should be equal to 1 to con- 
sider these.) 


Indicates whether there is a constant heat flux boundary 
condition along the larger tube. 


Not constant heat flux; the program needs a specified 
wall temperature distribution 


Constant heat flux boundary condition 


QADD 


TIN 


SWNO 


=0 
Zany value 


UFUIN 


RLARGE 
& 
DNOZLE 


IDATA 


A-3 


Wall heat flux in the case of a constant heat flux 
boundary condition, J/m@s. (To be used only if 
ICHF = 1) 


Inlet temperature, K° 


Swirl number, defined as the axial flow of tangential 


momentum divided by the axial momentum rate, 


SRitir?ar/e [¥i* rar, where R is the nozzle radius. 
No swirl 

Magnitude of swirl number 

The maximum inlet velocity (on the centerline); 
important only if SWNO = 0. The program then 


treats the inlet as fully developed pipe flow and 
UFUIN is the value at the centerline. 


DNOZLE /2 RLARGE 
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Parameters specify if the input velocities, turbu- 
lence kinetic energy and dissipation will be sup- 
plied via READ statements or will be calculated 

by the program as fully developed pipe flow. 


For swirl number = 0 (SWNO = 0) 


For others. In this case, at the University of 
Arizona Computer Center, the program Teach.Vel 
should be run first; its output should be saved 
on Tape 4; then the present_program is run and 
it will automatically read U, W, K, e from Tape 4. 


Number of nodes in x and y directions, respectively 


CRE Et Le 


Number of radial nodes in inlet tube. 
//////// 


JNOZLE——————> 


Locations of nodes in the x-direction 


Location of nodes in the radial direction (starts 
from centerline) 
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TWALN North wall temperature, K 


TWALW Side wall temperature, K 
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